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Recent observations by the CoGeNT collaboration (as well as long standing observations by 
DAMA/LIBRA) suggest the presence of a ~ 5-10 GeV dark matter particle with a somewhat large 
elastic scattering cross section with nucleons (a ~ 10~ 40 cm 2 ). Within the context of the minimal 
supersymmetric standard model (MSSM), neutralinos in this mass range are not able to possess 
such large cross sections, and would be overproduced in the early universe. Simple extensions of 
the MSSM, however, can easily accommodate these observations. In particular, the extension of the 
MSSM by a chiral singlet superfield allows for the possibility that the dark matter is made up of 
a light singlino that interacts with nucleons largely through the exchange of a fairly light (~30-70 
GeV) singlet-like scalar higgs, h 1 . Such a scenario is consistent with all current collider constraints 
and can generate the signals reported by CoGeNT and DAMA/LIBRA. Furthermore, there is a 
generic limit of the extended model in which there is a singlet-like pseudoscalar higgs, oi, with 
m ai ~ mhi and in which the x°X° an d bb, ss coupling magnitudes of the hi and a\ are very similar. 
In this case, the thermal relic abundance is automatically consistent with the measured density of 
dark matter if m x a is sufficiently small that X°X° bb is forbidden. 

PACS numbers: 95.35.+d, 12.60.Jv, 14.80.Da; FERMILAB-PUB-10-316-A 



Recently, the CoGeNT collaboration has reported the 
detection of very low energy events which cannot be ac- 
counted for with known backgrounds [l|. It has been 
shown than it is possible to interpret these events as 
the elastic scattering of a light dark matter particle 
(m ~ 5 — 10 GeV) with a cross section on the order 
of ~10~ 40 cm 2 pHa]- Intriguingly, the range of masses 
and cross sections implied by CoGeNT is not very far 
from the region required to explain the annual modu- 
lation observed by the DAMA/LIBRA collaboration @- 
8j . While null results from liquid XENON- based experi- 
ments Q and CDMS [l(| somewhat constrain dark mat- 
ter interpretations of the CoGeNT and DAMA/LIBRA 
signals, uncertainties in the dark matter velocity distri- 
bution, as well as the various experiments' energy scale 
calibrations and quenching factors (and/or scintillation 
efficiencies) could potentially lead to a consistent inter- 
pretation [BL ITij . 

Since the announcement of the CoGeNT result, a num- 
ber of groups have begun to explore the dark mat- 
ter phenomenology of this signal [!, [l2|-Q- Within 
the context of the Minimal Supersymmetric Standard 
Model (MSSM), dark matter explanations for the Co- 
GcNT/DAMA signals face considerable challenges. The 
range of elastic scattering cross sections predicted for 
neutralinos falls more than an order of magnitude short, 
even in the most optimal regions of parameter space fl5l . 
ITrjj . While this could plausibly be reconciled by adopting 
a significantly higher local density of dark matter (the re- 
quired cross section scales inversely with the local dark 
matter density), the relic abundance of very light (5-10 
GeV) neutralinos in the MSSM is also predicted to be 
well above the measured cosmological dark matter den- 



sity [la ] (for earlier work on light neutralino dark matter 
in the MSSM, see Ref. [13 )• Thus, even in optimistic 
regions of the MSSM parameter space, it is very diffi- 
cult to accommodate the observations of CoGeNT and 
DAMA/LIBRA. 

To increase the elastic scattering cross section and re- 
duce the thermal relic abundance of neutralino dark mat- 
ter, one could consider a combination of larger couplings 
or lower masses for the particles exchanged. However, for 
very light neutralinos (whose scattering is typically domi- 
nated by scalar higgs bosons) this prospect is constrained 
by Tevatron and LEP II data which require the masses 
of MSSM higgs bosons to lie above > 90 GeV. However, 
the constraints need not apply in supersymmetric mod- 
els with extended higgs sectors [18j . As we will show, 
in such scenarios it is possible for a 5-10 GeV neutralino 
to produce the observed signal through the exchange of 
a light (^30-70 GeV) scalar higgs, while also generating 
the correct thermal relic abundance. 

Generically, the spin-independent elastic scattering 
cross section of dark matter with a nucleus is written: 

a „ 4m 2 , m 2 , + _ 

tt(todm + m N y 

where tun is the mass of the target nucleus (of atomic 
number Z and mass A), and 771dm is the dark matter 
mass. f p and /„ are the dark matter's couplings to pro- 
tons and neutrons: 
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where a q are the dark matter's couplings to quarks (in 
the Lagrangian) and f^' n \ f^c^ are hadronic matrix 
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elements [lj|. An appropriate nuclear form factor ac- 
counts for the effects of finite momentum transfer. We 
use the values f^j = 0.020 ± 0.004, f$ = 0.026 ± 0.005, 



= 0.118±0.062, and = 0.84 |20|. The uncertain- 
ties in determination of hadronic matrix elements, which 
are due to the uncertainty in determination of ovjv , might 
reach approximately a factor of two and therefore result 
in a factor of two uncertainty in the dark matter cou- 
plings. 

For light MSSM ncutralinos, the neutralino-quark cou- 
pling is dominated by scalar higgs exchange (contribu- 
tions from squark exchange are typically negligible). For 
down-type quarks, this coupling is (20| : 
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where the Nu's denote the composition of the lightest 
ncutralino (x° = NuB + N r2 W a + Ni 3 H d + Ni 4 H u ), and 
s a and c a denote the sine and cosine of a, which relate 
the scalar mass and gauge eigenstatcs. The correspond- 
ing expression for up-type quarks is found by replacing 
cos/3 <-> sin/3 and Nu N13. 

The largest elastic scattering cross sections in the 
MSSM arise in the case of large tan /3 and sin(/3 — a) ~ 1 , 
significant N13, and relatively light m#. In this limit, 
the lighter higgs, h, is approximately standard model-like 
and the heavier H is approximately H d , and one finds 
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where the reference values for the Higgs mass, tan/3, N 2 3 
and c a have been chosen in the most optimistic way as 
discussed below. 

The higgsino content of the lightest neutralino is con- 
strained by the invisible width of the Z as measured 
at LEP, r?^ p = 499 ± 1.5 MeV. In contrast, the stan- 
dard model prediction for this quantity is slightly (1.4cr) 
higher, if^ = 501.3 ± 0.6 MeV Combining the 

measured and predicted values, we find a 2a upper 
limit of T z _> x o x o < 1.9 MeV. As T z ^, x o x a scales with 
[7Vf 3 — N 2 4 ] 2 , we can translate this result to a limit of 
\Nf 3 — N 2 4 \ < 0.103. For moderately large values of 
tan (3, the two higgsino terms do not efficiently cancel, 
requiring |7V 13 | 2 < 0.103. 

rriH and tan /3 are constrained by a number of mea- 
surements, including those of the rare decays t — > bH + , 
B s — > B — > tv, and direct limits on higgs pro- 

duction followed by A/H — s- r + r _ . While these lim- 
its vary somewhat depending on the precise values of 



the MSSM parameters adopted, in general they imply 
tan/3 < 20 - 30 for m H ,m A - 90 - 150 GeV (the 
strongest limits coming from the latest LHC results). 
Constraints from LEP II further require ran, a ^ 90 GeV. 
When these limits are taken into account, we find that 
cr x°p-n ~ 10 -41 cm 2 [TH, [l6j|, which falls short of that 
implied by the CoGeNT and DAMA/LIBRA signal by 
about an order of magnitude. 

Furthermore, light neutralinos in the MSSM are in- 
evitably predicted to freeze out with a thermal relic abun- 
dance in excess of the measured dark matter density. 
Indeed, there is a kind of inverse relation whereby for 
nriH ~ ttia near 90 GeV, appropriate Q x oh 2 from annihi- 
lation via the A to bb and t + t~ is roughly proportional 
to the inverse of <7 x o p n and only falls to a value of or- 
der Q x oh 2 ~ 0.1 when o x o vn is of order the maximal 
value indicated above. Given that current experimental 
constraints in the MSSM context do not allow one to 
achieve this maximal value, the MSSM inevitably leads 
to too large a value for f2 x o/i 2 . Increasing or mn 
worsens the situation. 

To increase the cross section beyond the range allowed 
in the MSSM, an obvious direction is to consider models 
with lighter higgs bosons. As the cross section scales with 
the inverse of the fourth power of the exchanged higgs 
mass, even modest reductions in tuh could increase the 
cross section to the levels required. As an example of a 
framework in which light higgs bosons are possible, we 
extend the MSSM by a chiral singlet superfield S, con- 
taining two neutral scalars H$ and A$ and a Majorana 
fermion S. The theory is described by superpotential [22[ 



^fi s S 2 + nH u H d + \SH u H d + ^kS 3 , (6) 



(along with the MSSM Yukawa interactions) and soft La- 
grangian 



C soft - v%S + B^H u H d + ^m 2 s \S\ 2 + ^B S S 2 

+\A x SH u H d + \A K S 3 + H.c. (7) 



(along with the MSSM ^4-terms). Specific implementa- 
tions of such a singlet typically involve a subset of these 
terms. For example, in the Next-to-Minimal Supersym- 
metric Standard Model (NMSSM) @|, a Z 3 symmetry 
is imposed which only allows the terms involving A, k, 
A\ and A K . Here, we do not tie ourselves to this particu- 
lar model, but instead consider the full range of terms 
as described in Eqns. (j6|) and (J7|), which we refer to 
as the Extended Next-to-Minimal supersymmetric Stan- 
dard Model (ENMSSM). 

The tree-level ncutralino mass matrix in the 
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B, W 3 , H u , H d , S basis is 
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where u u and are the up- and down-type higgs vevs 
and s is the vev of the singlet higgs. A light (< 10 GeV) 
neutralino consistent with LEP II chargino searches must 
be either mostly B or S. The B docs not couple to a 
mostly singlet higgs, and that case is thus similar to the 
MSSM. From here on, we focus on the case where the 
lightest neutralino is mostly S. One might imagine that 
the strict NMSSM would allow sufficient flexibility. It 
has been shown [H] that for the bulk of the NMSSM, af- 
ter imposing LEP and B-physics constraints the lightest 
neutralino is always bino-like and elastic cross sections 
as large as required by CoGeNT and DAM A/LIBRA are 
not possible. Nonetheless they are "only" a factor of 
10 too small (whereas in Ref. |25| the largest cross sec- 
tion found was a factor of 100 too small). An exception 
to this conclusion was identified in Ref. (26[, in which a 
very light (^GcV) singlet scalar higgs is able to generate 
the very large elastic scattering cross section required by 
CoGeNT and DAMA/LIBRA within the context of the 
NMSSM. These scenarios require a considerable degree 
of fine-tuning of the parameters. As we show, in the EN- 
MSSM, it is possible to find less fine-tuned scenarios with 
large elastic cross section and correct relic density when 
the light higgs has mass > 30 GeV even when tan/3 is 
small enough to evade the most recent LHC limits on 
Higgs bosons with enhanced bb couplings proportional to 
tan /3. 

We proceed by engineering the lightest neutralino to 
be mostly S, together with a light higgs that is predom- 
inantly singlet. The lightest neutralino will naturally be 
predominantly singlino provided the quantity \2ks + /is\ 
is much smaller than |^ + As|, Mi and M 2 . For example, 
for k ~ 0.45, s - 2 GeV, ~ 0, A - 0.01, tan/3 - 15, 
large Mi, large M 2 , and /i ~180 GeV, we find that the 
lightest neutralino is singlino-like (iVi 5 = 0.99) with a 
mass of approximately 5 GeV. 

The conditions under which the lightest higgs, hi, is 
mostly singlet are somewhat more complicated. A simple 
limit which leads to desired phenomena can be obtained 
for small A. In the limit A — >• 0, but keeping \A\ of mod- 
erate size, the singlet decouples from the MSSM (which 
has standard higgses), and has a mass determined pri- 
marily by the XA\v 2 sin2/3/(2s) singlet-singlet entry in 
the mass-squared matrix, with weak dependence on Bs, 
m|, us, k, and A K . Then, provided the soft terms for the 
singlet are sufficiently small, and Aw is much smaller than 
to/j, the singlet represents a perturbation on MSSM higgs 
phenomenology, with a light singlet state mixed with the 
MSSM to a degree controlled by A. In this limit, the light 
CP odd higgs will also be predominantly singlet, with 



mass and mixings to the MSSM A that are proportional 
to \A\v 2 sin2/3/(2s) with the same coefficients as for the 
hi. Effectively, the hi and ai combine to form a sin- 
gle complex nearly singlet scalar state. Through mixing 
with the MSSM higgses, both the hi and ai have cou- 
plings proportional to the usual MSSM interactions, but 
reduced by the small amount of mixing. For the param- 
eters listed earlier, along with A K ~ 33 GeV, A\ ~ 3400 
GeV, B s ~ 0, m| ~ one finds m hl ~ 40 GeV, 
m ai ~ 35 GeV with \F s (hi)\ 2 ~ \F s (ai)\ 2 ~ 0.65 and 
\F d (hi)\ 2 ~ \F d (ai)\ 2 ~ 0.35, where F s (hi) and F d (hi) 
are the singlet and components (at the amplitude 
level) of the hi and F 8 (ai) and F d (ai) are the singlet 
and A d components of the oi. To reemphasize, in the 
scenarios we consider the ai and hi are close in mass 
and have F s and Fd components that are very similar. 

The singlino coupling to down-type quarks via hi ex- 
change at Lagrangian level is given by 

Od_ = g 2 KN 2 5 F s (hi)F d (hi) 
m d 2y/2m w mf ii cos (3 

For hi exchange only, the resulting spin-independent 
cross section is 

a x o p , n « 3 .3xl0- 40 cm 2 ^ 

which is of order the value required by CoGeNT and 
DAMA/LIBRA. Furthermore, the mostly singlet nature 
(\F s (hi)\ 2 ~ 0.65) of the hi easily allows it to evade the 
constraints from LEP II and the Tevatron, as we dis- 
cuss below. The fact that moderately large \F d (hi)\ 2 is 
required argues that this scenario will be difficult to ar- 
range in the (MSSM) decoupling regime of m a . 2 3> mh 2 , 
implying that all of the mostly MSSM higgses are likely 
to have masses only slightly above the LEP II limit. In 
fact, for the parameters listed above, m,h 2 ~ 109 GeV 
(but escapes LEP limits since it has substantial singlet 
component, \F s (h 2 )\ 2 ~ 0.28), m a2 ~ 111 GeV and 
to/i 3 ~ m h + ~ 125 GeV. After including h 2 and /13 in the 
computation of o~ x o pn the coefficient of 3.3 is reduced to 
about 2.2 in Eq. (fTU|) due to some partial cancellation at 
the amplitude level. 

Collider constraints on the hi and h 2 are largely evaded 
for sufficiently small \F s (hi)\ and \F s (h 2 )\. LEP II places 
constraints through production of Zhi (for the scalars) 
and pair production of aihj for the pseudoscalars. The 
heavier mass eigenstates (/13, and a 2 ) look like their 
MSSM counterparts, with their couplings slightly re- 
duced by a small singlet component. Thus, provided they 
represent a viable point of MSSM parameter space, they 
will be allowed here as well. The light (mostly singlet) hi 
and ai must have small enough Z-Z-hi and Z-hi-ai inter- 
actions to be consistent with existing searches. In the sce- 
narios considered here, the coupling of the hi to WW, ZZ 
(relative to the SM coupling), denoted Cv{hi), is very 
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small, |Cy(/ii)| < 0.1, which easily allows ^40 — 50 
GeV to be consistent with LEP limits on Zhi [13 ■ Small 
|Cy(^-i)| arises in the limit A\ > /i. The Z-h\-a\ cou- 
pling is similarly suppressed implying that LEP limits on 
h\a\ pair production are easily evaded. Pair production 
of a\ together with the mostly SM-like light higgs /12 is 
sufficiently suppressed by the largely singlet nature of the 
01. The Tevatron can produce scalars and pseudoscalars 
through the reaction bg — > b + h, a, where the h or a can 
decay into either bb or t + t~ pairs. In our scenarios, the 
strongest constraints arise for h = h 3 and a = a 2 since 
they are fairly light and mainly non-singlet and have bb 
couplings that are enhanced at large tan/3. Null LHC 
searches with L = 1 fb _1 of accumulated luminosity at 
CMS require tan (3 < 25 for a doublet-like h^oi 02 with 
(non-degenerate) masses of order 100 GeV 28j. A sim- 
ilar limit on tan/3 can be obtained from the null search 
for t —} H + b [23]. Taken all together, the central pa- 
rameters of Eqn. (|10j) are on the border of a number of 
higgs searches, and are thus being tested by end-phase 
Tevatron and early LHC running. 

The thermal relic density of neutralinos is determined 
by the annihilation cross section and neutralino mass. 
In the m x o range we are considering here, the poten- 
tially important annihilation channels are to bb or t + t~ 
through the s-channcl exchange of Higgs bosons. Given 
that m/jj ~ m ai and that the F s and components of 
the h\ and a\ are similar, the CP-odd a\ is dominant, 
annihilation via the CP-even h\ being p-wave (v 2 ) sup- 
pressed. 1 The annihilation cross section for the bb final 
state that results from s-channel exchange of the a\ is 
given by 

7V c52 2 K 2 m 2| Fs(ai) |2| Fd(ai) |2 



47rm 2 1/ cos 2 (3 

to 2 (1 - ml/m 2 x0 ) 1 / 2 

X (4m 2 - m 2 ai ) 2 + m 2 ai T 2 ai 



(11) 



where v is relative velocity between the annihilating neu- 
tralinos, N c = 3 is a color factor and r ai is the width of 
the exchanged higgs. The annihilation cross section into 
t + t~ is obtained by replacing mj, — > m T and iV c — > 1. 
The thermal relic abundance of neutralinos is obtained 



as 



Mpi Tpoy/gl {a x o x ov} 



(12) 



where g+ is the number of rclativistic degrees of freedom 
available at freeze-out, (a x o x ov) is the thermally aver- 
aged annihilation cross section at freeze-out, and Tfo is 
the temperature at which freeze-out occurs. In the sce- 
narios considered here, one finds that the annihilation 
rate is too large if the a\ — > bb channel is open. As a 



1 This differs from the scenario of in which the a\ is highly 
singlet and annihilation is dominated by the h\. 



result, consistency with the observed relic density is only 
"automatically" obtained if m x o < m^pole) ~ 5.28 GeV. 

For the range of masses and cross sections considered 
here, we find m x o /T-po ~ 20, leading to a thermal relic 
abundance from a\ — > t + t~ of 
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applicable so long as m x o is small enough that the bb 
channel is not open. After including other exchanges, 
one obtains a value of Q x oh 2 that is consistent with 
the measured dark matter density, f^cDM^ 2 = 0.1131 ± 
0.0042 H|. 

Thus, the desired relic density automatically results 
at low m x o once the relevant combination of couplings 
and higgs masses are set to accommodate CoGcNT and 
DAMA/LIBRA. This confluence of parameter space is 
peculiar to models with scalar/pseudoscalar exchange 0, 
l3ll [32( 1 . For example, a Dirac fermion or a scalar with 
vector interactions will cither overproduce the CoGcNT 
and DAMA/LIBRA rates or will predict a thermal relic 
density in excess of the measured dark matter abundance. 

This scenario also has interesting implications for the 
indirect detection of dark matter. In particular, as the 
dark matter annihilation rate in any given region scales 
with the inverse of the square of the dark matter mass, 
the light neutralino we are considering could, in princi- 
ple, lead to enhanced fluxes of various annihilation prod- 
ucts [33l ]. Quantitatively, the spectrum of gamma-rays 
from dark matter annihilations can be written as 



$ 7 ( J B 7 ,V) 
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p 2 (r)dl, 



(14) 



where av is the dark matter annihilation cross section 
multiplied by the relative velocity of the two neutrali- 
nos, ip is the angle observed relative to the direction of 
the Galactic Center, p(r) is the dark matter density as a 
function of distance to the Galactic Center, dN 1 /dE 1 is 
the gamma ray spectrum generated per annihilation, and 
the integral is performed over the linc-of-sight. For a neu- 
tralino with m x o ~ 5 GeV and with an annihilation cross 
section of av as 3 x 10~ 26 cm 3 /s to t + t~, the annihila- 
tion rate in the Galactic Center (assuming a NFW halo 
distribution) is predicted to lead to a flux of gamma-rays 
above 1 GeV of as 2.9 cm -2 yr _1 from the inner degree 
of our galaxy, corresponding to thousands of events per 
year observed by the Fermi Gamma Ray Space Telescope 
(FGST). Furthermore, the flux and spectral shape of the 
gamma-ray emission observed by the FGST from this re- 
gion of the sky is quite similar to that predicted from dark 
matter annihilations [351. FGST's observations of dwarf 
spheroidal galaxies [34| are also potentially sensitive to a 
dark matter particle with these characteristics. 

The prediction that the dark matter annihilates pri- 
marily to t + t - also insures that our dark matter can- 
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didate will not violate constraints from cosmic ray an- 
tiproton measurements, as set by the PAMELA experi- 
ment [13, HE • If gamma ray searches for the prod- 
ucts of dark matter annihilations were in the future to 
constrain the low-velocity annihilation cross section to 
be well below the value of av ~ 3 x 10~ 26 cm 3 /s, this 
would rule out the present model as well as the possibil- 
ity that the dark matter is a scalar with scalar interac- 
tions urnm. 

Another indirect detection mode with decent prospects 
is to search for energetic neutrinos produced through the 
capture and annihilation of dark matter in the core of 
the Sun. The Sun is predicted to capture dark matter 
particles at a rate given by 



tion of their rate at the level of 2.8a [4l|. This result 
provides further motivation for the type of model consid- 
ered in this paper. 
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where p x o is the local dark-matter density, v is the local 
root-mean-square velocity of halo dark-matter particles, 
and cth and one are the elastic scattering cross sections of 
the WIMP with hydrogen and helium nuclei, respectively. 
In the model under consideration, the elastic scattering 
cross section of the neutralino is sufficiently large that 
the processes of capture and annihilation quickly reach 
equilibrium in the Sun, removing any dependence on the 
neutralino's annihilation cross section. The high capture 
rate in this model is predicted to produce a sizable flux 
of GeV-scale neutrinos, comparable to the constraints 
currently placed by Super-Kamiokande 0, [H, H(j. In 
particular, for m x o ~ 5 GcV and annihilation entirely to 
t + t~, Super-Kamiokande data can be used to constrain 
(J x°p ~ ^ x 10~ 41 cm 2 Q, for reasonable astrophysical 
assumptions. Larger volume neutrino experiments such 
as IceCube have energy thresholds which are too high 
to observe the annihilation products of such light dark 
matter particles. 

In summary, we have considered the possibility that 
neutralino dark matter is responsible for the signals re- 
ported by the CoGcNT and DAMA/LIBRA collabora- 
tions. Although, the elastic scattering cross section of 
neutralinos with nuclei in the MSSM is too small to ac- 
count for these observations, the same conclusion is not 
necessarily reached in extended supersymmetric models. 
In particular, we have discussed models in which the 
MSSM is extended by a chiral singlet supcrfield. In such 
a model, a light singlino-likc neutralino, which interacts 
with nuclei through the exchange of a largely singlet- 
like, scalar higgs, can possess an elastic scattering cross 
section capable of generating the observations reported 
by CoGeNT and DAMA/LIBRA. Furthermore, the sce- 
narios considered automatically lead to a thermal relic 
abundance of neutralinos consistent with the observed 
density of dark matter for m x o < 5 GeV. 

After the completion of this project, the CoGeNT col- 
laboration reported the detection of an annual modula- 
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